Chain deformation helps translocation 
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Deformation of single stranded DNA in translocation process before reaching the pore is investigated. By solving the Laplace 
equation in a suitable coordinate system and with appropriate boundary conditions, an approximate solution for the electric field 
inside and outside of a narrow pore is obtained. With an analysis based on "electrohydrodynamic equivalence" we determine the 
possibility of extension of a charged polymer due to the presence of an electric field gradient in the vicinity of the pore entrance. 
With a multi-scale hybrid simulation (LB-MD), it is shown that an effective deformation before reaching the pore occurs which 
facilitates the process of finding the entrance for the end monomers. We also highlight the role of long range hydrodynamic 
C*~) interactions via comparison of the LB-MD results with those obtained using a Langevin thermostat instead of the LB solver. 
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1 Introduction 

The detailed dynamics and physics of polymer translocation 
through confined geometries, e.g. channels and pores, are not 
well-understood 1 ' 2 . In the last decades, experiments with sin- 
gle nanopores have been used to study translocation dynamics 
of single- and double-stranded DNA macromolecules through 
pores 3 . In these studies, a voltage bias drives the macroions 
through a single nanopore which can be a protein channel or a 
pore on an artificial solid-state substrate 4-6 . When a polymer 
traverses a pore in a membrane, it partially blocks the chan- 
nel and so decreases the magnitude of the ion current passing 
through the pore due to an applied voltage. By recording the 
ionic current modulation, one can thus monitor the polymer 
translocation indirectly 1 ' 3 ' 5 ' 7 . 

Due to the tiny size of the pore in comparison with the elec- 
trolyte chamber, the electrostatic voltage drops mainly inside 
the pore. Based on this observation, many studies neglect the 
effect of the electric field outside the pore on the dynamics of 
translocation 8 " 10 . In such studies, it is also often assumed that 
polymer reaches the pore by a purely thermal diffusion process. 
Once sufficiently close to the pore, one of the end monomers 
then has to find the pore entrance. This marks the beginning 
of the translocation process during which the polymer is drawn 
into the pore by the electric force. However, the electric field 
in the vicinity of the pore entrance may have an important ef- 
fect on DNA capture. In some studies in which the capture rate 
of translocation phenomenon is investigated, the biased diffu- 
sion of polyelectrolyte in a region close to the pore is consid- 
ered 2 ' 311 " 13 . In this paper, we show that this is not the only 
effect of the extended electric field outside the pore on translo- 
cation. Electric field is highly convergent on the cis side of the 
pore and divergent on the other side 2 ' 12 . The fact that DNA is a 
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highly charged flexible macromolecule, is a strong incentive to 
study the possibility of DNA deformation and polarization in a 
nonuniform electric field, which is prevalent close to the pore 
entrance. There is experimental evidence that, before translo- 
cation, DNA deformation and counterion decoupling occur in 
a capture region comparable to the hydrodynamic size of the 
chain 14 . It is argued that, upon increasing the voltage, this 
mechanism can enhance the rate of diffusion-limited capture 
of DNA 1 1 . Deformation of DNA and its effect on the capture 
rate of the translocation process has been addressed in a num- 
ber of previous theoretical works 12 ' 1315 . To the best of our 
knowledge, however, this issue has not been tested in a direct 
simulation study so far. 

Applying a mechanical 1617 , electrical 18-20 or hydro- 
dynamical 21 " 23 force on a long polymer (in comparison to its 
persistence length) deforms the relaxed configuration of the 
polymer 24 ' 25 . According to this effect, a polyelectrolyte can 
be deformed from its coiled equilibrium configuration to an 
extended chain by a uniform (for an anchored polymer) or a 
non-uniform (for a freely moving polymer) electric field 26 ' 27 . 
By simulating the DNA motion in a network of insulated cylin- 
ders, Tessier and coworkers have shown that DNA can deform 
electrophoretically in the entropic traps 28 . Randall et al. have 
also observed free DNA deformation in an electric field gradi- 
ent. These authors have shown that the electric field gradient 
created by insulator obstacles deforms DNA molecules. Us- 
ing Long's electrohydrodynamic equivalence 29 , they have also 
interpreted the observed dynamics of DNA deformation in an 
electric field gradient via a kinetic model 26 . 

Recently, Kowalczyk et al. used an oblate spheroidal co- 
ordinate system to solve the Laplace equation inside the pore 
and have determined the electric field as a function of the pore 
shape parameters and the applied voltage 30 . This assumption 
yields a reasonable description of access resistance of a pore 
before and after the insertion of a DNA molecule into the pore. 
In this paper, we use this coordinate system and, after a sim- 
plification of the geometry of the pore, find an approximate ex- 
pression for the electric potential in the entire space, just by 
setting the electrodes at infinity. We show that the correspond- 
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ing electric field is consistent with recent experiments on the 
conductivity of nanopores 30 . We then determine the velocity 
gradient tensor and its eigenvectors and discuss the possibility 
of DNA deformation due to the presence of electric field gra- 
dient near the pore entrance. Such deformation can help the 
ends of DNA molecules find the pore entrance and thus facil- 
itate the capture process. Indeed, experiments show that DNA 
is more likely to be captured by the pore at the DNA's ends 31 . 
A very recent study confirms these findings by a more accurate 
statistical analysis of the data 32 . 

Regarding the simulation methodology, we employ a hybrid 
lattice Boltzmann - molecular dynamics (LB-MD) technique. 
The coupling to LB is motivated by the need to account for the 
effect of long range (~ 1/r) hydrodynamic interactions (HI). 

For a single-stranded DNA (ssDNA) in dilute solutions 
which is studied in this paper, the Debye length, Ad, and per- 
sistence length, / p are of the same order of magnitude (for ex- 
ample, a typical value of Ad in a 10 mM solvent of monovalent 
salt is 3 nm or larger while / p ~ lnm 27,33 ). So the local elec- 
tric field disturbances by DNA and counterions movement can 
not be ignored 33 . On the other hand, it seems that, even if we 
ignore this effect, we cannot neglect the importance of the ar- 
rangement of counterions around the DNA chain when study- 
ing the dynamics of a poly electrolyte in a non-uniform electric 
field. The ability of counterions to glide along the chain and 
detach and reattach to it, can give rise to a local polarization 
of the chain and thereby induce an extra interaction between 
the charged polyelectrolyte and the external field. Therefore, 
we also study the effect of counterion rearrangement in a inho- 
mogeneous electric field in the simulation part of this paper by 
applying a P3M method. 

For simplicity, we neglect the effect of image charges and 
the related dielectric mismatch between the membrane and the 
solution. This is a good approximation if the Debye length 
is below the shortest characteristic distance between the DNA 
molecule and the membrane 33 . 



2 Theory and Methods 
2.1 Electric field 

To find an approximation for the electric field of a pore in 
the whole space, we use the oblate spheroidal coordinate sys- 
tem 30 but with the electrodes at infinity. We use ju, v and (j) 
as the oblate spheroidal coordinates with a focal ring of radius 
c, where fi e (-oo 7 +oo), v g [o,?r] and e [0,2n) (Fig. la)). 
These selected coordinates map to the cylindrical ones by the 
simple relations: p = ccosh/icos v, z = csinhjtisinv and <j> = 
0. Surfaces of constant Li are confocal semi-ellipsoids, sur- 
faces of constant v are one-sheeted hyperboloids of revolution 
and surfaces of constant (j) are semi-planes normal to the x-y 




Figure 1 a) One can mathematically simplify the geometry of the pore 
by considering it as a surface of a hyperboloid with v = V(>. b) The 
DNA translocation process is shown schematically, in the presence of 
an extended electric field (blue lines). DNA monomers are represented 
in black and some of its counterions in red. 



plane. This set of coordinates map to the Kowalczyk's notation 
with ji = sinh -1 1, v = sin -1 rj; is the same in the two nota- 
tions. We assume that when the pore radius is greater than its 
length (a > I, where a and / are the pore radius and length, re- 
spectively), we can substitute a one-sheeted hyperboloid with 
V = Vo for the membrane wall and the pore (Fig. la), where 
Vo can assume any arbitrarily small value (Vo <C 7t/2) which 
can be obtained from the pore radius, a, by the simple relation, 
Vo = cos - (a/c). Since we want to study the deformation of 
a negatively charged polymer, it is preferable to consider the 
direction of electric field lines from the anode to the cathode 
(Fig. lb). We adhere to this convention throughout this paper. 
Assuming the cathode and anode sufficiently far from the pore 
center (negative and positive infinity, respectively), we have the 
following two boundary conditions: 



<&(jU,V,0)| M = ±O e = T 



(1) 



E(n,v,(j)).n\ v=Vo = V<J>(ju,v,0)| 



v=v 



0. 



Solving the Laplace equation with these boundary conditions, 
one obtains the following continuous electric potential in the 
entire space: 



Vo 



<£(/i,V,0) = — tan '(sinhju) 



(2) 



This potential is almost constant far away from the pore but 
varies rapidly close to it. The equipotential surfaces are con- 
focal semi-ellipsoids of constant Li with azimuthal symmetry. 
The field lines lie on the hyperboloids, confocal with the hy- 
perboloid representing the membrane wall (Fig. lb), 
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Ka cosh ji \J sinh 2 \i + sin 2 v 
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where is the coordinate scale factor in fx direction 
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Although we have assumed a > I to justify our simplification 
of the pore geometry, it is worthwhile to mention that one can 
construct pores with inverse aspect ratio (Z > a) by adding two 
one-sheeted semi-hyperboloids (with Vo <C 1) to the end of a 
cylinder of optional length (7 C ). Then, according to the unique- 
ness theorem, the solution of the Laplace equation is a continu- 
ous piecewise differentiable function comprised of three parts: 
first in the region z > l c /2, the electric field can be obtained 
from Eq. (3) after a translation by an amount of l c /2 in z di- 
rection; for —l c /2 < z < lc/2, the electric field is constant and 
homogeneous; finally for z < —k/2, the electric field can again 
be obtained from Eq. (3). We shall be mainly concerned with 
the chain dynamics before it reaches the pore. Thus, the details 
of the pore geometry are expected to have a minor effect on 
our results. We will return to a discussion of this aspect in the 
following sections. 

Near the pore entrance, the functional form of the electric 
field is in good agreement with Nakane's solution 2 . For r^§> a, 
the electric field lines approximately lie on the lines asymptotic 
to the hyperboloids which are parallel to the f direction. Hence, 
far away from the pore, we reach the radial approximation for 
the electric field 1 1 . 

Using this extended electric field, we can calculate the con- 
ductance of the pore by computing the current obtained from 
the integration of the current density, j = a s E, on any surface 
crossing the pore (a s is conductivity of the solution). Equipo- 
tential surfaces with constant jj, are good choices for this pur- 
pose. However, the special case of the jj, —0 surface turns out 
to be still more convenient. Using Eq. (2) and after some sim- 
plification, one obtains: 



J J-P-ds 



k/2 f 2n d ® u h 



vo 



d[i hix 



= 2g s Vqo 



1 — sin Vo 
cos Vo 



(4) 



From this equation one can calculate the conductance of the 
semi-space from infinity to the entrance of the pore, G acC ess = 
4(7 s fl(l — sinvo)/(cos Vo). As discussed above, for the general 
case one can consider the pore as a combination of a cylinder 
and two semi-spaces with membrane walls coinciding with a 
hyperboloid. Conductance of such a system is 



1 + 

^channel G"a 
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2 A a cos Vo x -i ... 

O s a z (- + -- ) . 5) 

v % 2 1 - sin Vn ' 



For Vo <C Tl/2 this equation approaches the Hall approximation 
for the pore conductance 34 , which provides an upper limit for 
the access resistance 30 . It has been shown in previous studies 
that the asymptotic behavior of G for the both limits of a <C / 
and a ^> / is in good agreement with experimental measure- 



2.2 Velocity gradient tensor 

To describe DNA deformation outside the pore we adopt the 
method introduced by Randall 26 . In our case, we generalize 
the method to three dimensions 35 . Along this line, we first 
note that the symmetric rate-of-strain tensor governs the local 
deformation of fluid elements solely. The local deformation in 
3D extensional flows can be represented by a 3 x 3 symmetric 
matrix with orthogonal eigenvectors and real eigenvalues. An 
ideal flexible object will exponentially extend/compress along 
the axis of extension/compression (with eigenvectors corre- 
sponding to positive/negative eigenvalues) at rates equal to the 
eigenvalues of the deformation tensor. This ideal deformation 
is known as affine deformation 26 . The symmetric matrix of 
deformation is equal to the symmetric part of the velocity gra- 
dient tensor. Electrophoretic velocity of a charged particle with 
mobility coefficient /i e i in an electric field E(jU,V,0) is equal 
to jj. e \E(jj., V, 0). Since V x E = 0, the gradient of this velocity 
vector field is symmetric, so its eigenvalues directly give us the 
strain-rate tensor. The velocity gradient tensor can be written 
as: 

/ABO 
jU e iV£ = S [ B C 
\ D 

where according to Eq. (3) we have, 



(6) 



a 2 n (sinh 2 jx + sin 2 v) 2 cosh 2 jx' 
A = +sinhjU(cosh2jU + sin 2 v), 
B = +cos v sin v cosh ju, 



(7) 



C = — sinh jU cosh fi, 

D = — sinh jU (sinh 2 ju + sin 2 v) . 

Now we can obtain the electrophoretic strain rates e (r) by cal- 
culating the eigenvalues 



. A+C±^(A-C) 2 +4B 2 
A± - S , 

A* = SD, 
and the corresponding eigenvectors, 

e± = (^-C)A+Bv, 



(8) 



(9) 
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of the tensor in Eq. (7). The eigenvectors specify the extension 
and compression axes of deformation. e± are normal to and 
lie in the jx-v plane. For test purpose and without any loss 
of generality, we restrict ourselves to the pore axis, defined by 
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Figure 2 a) Relative strain rate and relative electric field magnitude, 
normalized to their maximum values, versus relative distance to the 
pore, normalized to the pore radius (computed on the z axis. See 
Eq. (11)). b) Direction and relative magnitude of extension axis along 
the centerline of the pore. 

V = k/2 (z axis) 26 . On this axis, the eigenvalues are: 



H-Iv=jt/2 



= +2 



VjJ. e i sinh/z 
a 2 n cosh 2 ^ 



(«+=A), 



V=7T/2 



VjXgi sinh^ 



(10) 



V,e<4 



a 2 7t cosh 2 ^ 

As seen from Eq. (10), there is symmetry with respect to rota- 
tions around the pore axis. A similar symmetry is thus expected 
with regard to deformations in the plane perpendicular to the 
pore axis. On the other hand, a survey of the eigenvalues and 
the corresponding eigenvectors reveals that, in the cis cham- 
ber (sinh fi > 0), DNA elongation occurs along the fx direction 
while compression takes place in the normal plane. To sim- 
plify the equations, we write the extension eigenvalue (strain 
rate along the pore axis) in Cartesian coordinates: 



alt cr 



(11) 



A polymer along the centerline of the pore experiences a 
pseudo-homogeneous extensional field 26 . Fig. 2a shows the 
strain rates (normalized to the maximum value) and the mag- 
nitude of the electric field along the z axis. It can be seen that, 
while the maximum value of the electric field occurs exactly 
at the middle of the pore, the strain rate (and so the gradient of 
the electric field) has a maximum at a somewhat larger distance, 
which depends on the pore shape. The direction of the exten- 
sion axis and relative elongation along the centerline of the pore 
are also shown schematically in Fig. 2b. The above discussion 
suggests that the linear alignment of a polyelectrolyte chain is 
enhanced by the electric field. 
The integrated strain, 



e[x(t)]dt, 



(12) 



for a point-like molecule approaching the pore along the center- 
line can be calculated using Eq. (11). This enables us to com- 
pute e for a single chain using its center-of-mass trajectory, and 



draw the extension-strain curves and compare theoretical pre- 
dictions with the results of simulations of the next section 26 . 



2.3 Simulation 

The translocation of a ssDNA through a single nanopore, in 
the presence of an applied electric voltage, is simulated by em- 
ploying the ESPResSo package 36 . To take into account both 
the hydrodynamic and electrostatic interactions, we use a hy- 
brid mesoscale method in which the lattice Boltzmann (LB) 
method 37 " 42 is applied to simulate the hydrodynamic effects in 
the system and a P3M algorithm is used to calculate full elec- 
trostatic interactions between charged monomers and counte- 
rions 43 ' 44 . A bead-spring coarse-grained approach, in which 
groups of atoms are represented by a single interacting parti- 
cle, is used for the MD part of the simulation. This enables a 
computationally efficient study of the translocation process. In 
experimental setups for the investigation of translocation, salt is 
always added to the system. Although most of our simulations 
have been performed in the absence of salt, we have also sim- 
ulated a number of cases with added salt so that a conclusion 
can be drawn as to the effect of salt (see below). 

The pore is represented by a solid substrate containing a cir- 
cular orifice with symmetry axis along the z direction. The in- 
teractions between beads as well as between the chain and the 
membrane are written as the sum of five terms: 



U (r^) = t/FENE + t/wCA + tAvall + U e l ec + t/ e xt, 



(13) 



where I/fene and t/wcA ar e respectively, the chemical bonds 
and excluded volume potential between the chain beads. £/ wa n 
is a hard-core interaction between the chain monomers and 
membrane (including the pore). The last two terms represent 
the electrostatic interactions and the interactions with the ex- 
ternal electric field. These are written individually as 45 : 



t^FENE : 

Uy/CA : 

tAvall : 
t/elec : 



E ^ 2 ln(l 



bonds 

EE T 4e([( 



-) 2 ), 



' ' ail ) n -{^f] + \), 



°^) 12 -(— ) 6 ]+l), 

' iw y ixv ~r 



(14) 



where r,-; = 

and r, 
l 



- fj | is the distance between the sites i and j 
is the distance between the site i and the wall. a,j = 



k(di + dj) is the arithmetic mean of the radii of the beads in- 
volved in the interaction. The f superscript indicates that the 
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summation is restricted to interparticle distances below the cut- 
off radius (ry < r c and r,„. < r c ). We use the monomer diameter, 
a, as the unit of length. The unit of energy is chosen to be e. 
The time is given in units of Tlj = 0\Jmj e, where m denotes 
the mass of a monomer. Unless explicitly stated, all the quan- 
tities in this paper are given in these reduced units (r.u.). Ions 
(counterions and salt ions) have a radius of d[ on — 0.425 (7 33 . 
In the FENE potential, the bond stiffness and the maximum al- 
lowed bond length are set to k — 30 and R = 1.5, respectively. 
We use a cut-off distance of r c = 2 1//6 (7, ; - to consider only the 
repulsive part of the Lennard-Jones potential. Together with 
the FENE potential, this leads to a bond length of rt, nd = L2. 
All simulations are performed at a temperature of k%T = 1 (the 
Boltzmann constant is set to unity (k-Q = 1))- tw is the wa- 
ter permitivity. The external electric field, E(r), is defined by 
Eq. (3). In all the studied cases, we add an appropriate amount 
of monovalent counterions to the charged polyelectrolyte in or- 
der to neutralize the system. To reach a desired salt concentra- 
tion, we add a corresponding number of coions and counterions 
to the system 33 . A P3M method implemented in the ESPResSo 
package is used to calculate the electrostatic potential between 
charged beads. The screening effect of counterions is thus con- 
sidered explicitly in the model. 

An approximate map of the present model to a physical sys- 
tem leads to the real units of e = 0.59 Kcal/mol, a = 1.0 nm 
and Tlj = 2 x 10~ n sec. 

Choosing e (the elementary charge) as the electric charge 
unit, the real unit of voltage is ~ 25 mV. In our simulations, 
depending on the situation, the applied external voltage varies 
between 8 to 12. Using the above unit conversion, this cor- 
responds to 200-300 mV which lies nearly at the middle of 
the range in typical translocation experiments with solid-state 
nanopores. Each monomer carries a reduced electric charge of 
—3 and each counterion and salt ion has a charge of ±1. The 
Bjerrum length Iq = e 2 /4-neoe w ksT = 0.70 corresponds to the 
water Bjerrum length at room temperature 45 ' 46 . 

The size of the simulation cell L z = 5Rg x 5Rg x 

IQR S , where R g denotes the gyration radius of the chain. The 
pore has a radius of 1.5(7 and an axial extension of 5(7. The 
pore radius is sufficiently large to allow the chain beads to pass 
through and small enough to prevent simultaneous passing of 
more than one bead (hairpin translocation). 

Most of the simulations reported here include hydrodynamic 
interactions via a frictional coupling of the particle dynamics 
to the LB method. We use the D3Q19 model implemented in 
the package with a kinematic viscosity of v = 3.0 and a fluid 
density of p = 0.864. The LB grid spacing is «lb = 1 and the 
LB time step is Tlb = 0.05. The coupling constant is set to 
Tbare = 20. 33 ' 46 . Again, all these quantities are given in the 
reduced (LJ) units. 

Additionally, random fluctuations for particle and fluid act 
as a thermostat. We performed a number of benchmark sim- 
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Figure 3 Typical configurations of a chain are shown in different sit- 
uations. Upper panels (a and b) are selected from the snapshots of a 
successful translocation. Extension along the electric field lines can be 
observed. Lower panels (c and d) show exactly the same situation for 
an unsuccessful try. Length of the chain is 30 beads and the simula- 
tions have been done with a reduced voltage of V* = 10 and no added 
salt. 



ulations to make sure that this model with the mentioned pa- 
rameters can reproduce the expected behavior of a single chain 
with HI. The results are in good agreement with those studies 
which have investigated the ability of MD-LB hybrid method 
in detail 47 " 49 . 

The applied electric field is turned on after an initial equili- 
bration run of duration t = 2 x 10 6 MD steps (At = 0.01). This 
time is sufficiently large for the chain to reach its equilibrium 
configuration for all investigated chain lengths, N. The poly- 
mer's center of mass is initially located on the pore's axis at a 
distance equal to the capture radius, r* , from the pore entrance. 
This is defined as the distance at which the electrophoretic in- 
teractions start to dominate and the polymer is attracted toward 
the pore, r* depends on the parameteres of the system such as 
the applied voltage, temperature, pore radius and the polymer 
length. It is obtained via a previously introduced method 12 and 
is found to lie between a few up to several times R g 12 . 

When DNA's center of mass reaches a distance greater than 
r*, we transfer it radially toward the pore by applying a decay- 
ing radial force which acts on the monomers only for distances 
r > r* and is zero otherwise. 

The DNA deformation, before entering the pore, is studied 
for some selected values of the electric voltage for chain lengths 
of N = 10,20, ...,60 in the absence of salt. For comparison, we 
have also performed a simulation for a 50 mM salt solution at 
a chain length of N = 20 monomer units. 
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Figure 4 Extension-strain 
curves for different volt- 
ages during the transloca- 
tion of DNA into the pore 
along the centerline of the 
pore. Each curve repre- 
sents the trajectory of a 
single molecule. Curves 
switch from solid black to 
dashed gray when the mid- 
dle monomer of the chain 
crosses the midpoint of the 
pore. Affine extension 
(e e , solid line) and com- 
pression (e~ E , dashed line) 
scaling are shown for ref- 
erence. In these figures, 
extension x ex is scaled to 
the polymer L and strain is 



scaled to 



4tl7l 



At. The ex- 



ternal applied voltage in- 
creases in successive fig- 
ures: a) V* =8b) V* = 10 
and c) V* = 12. N= 20 
and there is no added salt. 



3 discussion 

Due to the electrophoretic forces acting on the charged 
monomers, the molecule is attracted toward the pore. At the 
same time, the chain is deformed as a result of the non-uniform 
nature of the electric field. Upper panels of Fig. 3 show rep- 
resentative configurations of a chain before and during a suc- 
cessful translocation respectively. As can be seen, the polymer 
experiences a significant deformation before reaching the pore 
in such a way that one of the end monomers is captured by the 
electric field of the pore making the translocation possible. 

Even in an unsuccessful try (lower panels of Fig. 3), de- 
formation occurs before reaching the pore and on average the 
chain elongates along the field lines. This increases the prob- 
ability for an end monomer to be close to the pore entrance 
so that finding the entrance becomes easier for the chain in a 
subsequent try. 

One can sketch the extension-strain curves by calculating the 
integrated strain (Eq. (12)) and the chain's longest extension in 
each time step as described in previous sections. Relying on the 
assumption of affine deformation 26 , we expect a linear relation 
between the logarithm of the chain extension and the integrated 
strain. Such behavior can be observed in Fig. 4. In this figure, 
external voltage increases from V* — 8 to V* = 12. We have 
followed the successful translocations trajectory for those DNA 




10 15 
r/c 

V<w»v . w =60 
N=50 
. N=40 

• N=30 

• N=20 
> N=10 



30 




Figure 5 a) The radial ex- 
tension of the chain versus 
radial distance to the pore 
(in units of a), for differ- 
ent lengths of the chain. 
V* = 10 and there is no 
added salt. Inset shows 
| cos(0)| (0 is the angle be- 
tween the extension axis 
and the electric field direc- 
tion in the position of the 
polymer's center of mass) 
versus the radial distance. 
Data are averaged over 
tiny windows of radial dis- 
tance and over 100 simu- 
lations for each length, b) 
The same as (a) but with 
the radial electric field ap- 
plied outside the pore. 



molecules, whose path are nearly along the centerline of the 
pore. 

As can be seen, the chain initially extends in the z direction 
before reaching the pore. On average, the maximum exten- 
sions occur just before translocation and coincide nearly with 
the maxima of extension-strain curves. Just after the transloca- 
tion where the sign of the eigenvalues of velocity gradient ten- 
sor changes, we expect a compression along the z axis. This 
compression continues until the chain relaxes to its equilib- 
rium configuration and then the strain rates and integrated strain 
magnitudes go to zero. 

To see the effect of this deformation on translocation phe- 
nomenon, we computed the eigenvectors and eigenvalues of the 
polymer gyration tensor during the motion of the polymer. The 
eigenvector corresponding to the largest eigenvalue of this ten- 
sor is the extension axis of the polymer. Radial extension is 
defined as the projection of the end-to-end distance vector on 
the radial vector which points from the pore to the center-of- 
mass of the chain. In Fig. 5a, the radial extension of the chain 
(averaged over 100 independent simulations), is plotted versus 
distance to the pore. One can see from this plot that the exten- 
sion of the chain increases as the chain approaches the pore. 

To investigate whether this extension is along the electric 
field lines, we compute the angle, 9, between the extension 
axis of the polymer and the electric field direction at the posi- 
tion of the polymer's center of mass. The inset in Fig. 5a shows 
| cos(0) I (averaged over 100 independent simulations) versus 
distance to the pore. One can see that the elongation along the 
field lines is nearly perfect in the vicinity of the pore for all 
chain lengths investigated. 

Interestingly, we observe that the details of the emerging 
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• With salt (50 mM) 




■ No added salt 



















Figure 6 The radial ex- 
tension of the chain ver- 
sus the radial distance to 
the pore (in units of cr), 
with and without added 
salt. N = 20andV* = 10. 




Figure 7 a) Monomer approaching probability (MAP) along the chain 
for different lengths with explicit counterions for a voltage of V* = 10 
and no added salt. For better comparison between different lengths, 
we normalized the monomer ID (m) to the degree of polymerization 
(N). b) Same diagrams as in a) but using mean-field effective charges 
instead of explicit counterions. 



electric field only play a minor role for the behavior of the chain 
outside the pore. To illustrate this aspect, we show in Fig. 5b 
results on the chain deformation for a inhomogeneous but ra- 
dial electric field 1 1 . As a comparison of the panels (a) and (b) 
in Fig. 5 reveals, the deformation behavior of the polymer is 
nearly the same for the both types of electric field investigates. 

Regarding the effect of added salt, we investigated the elon- 
gation behavior of a chain in a solution of mono-valent salt with 
a concentration of 50mM (N = 20) for the two situations with 
and without added salt. As shown in Fig. 6, the extension of 
the chain along the electric field does not seem to be strongly 
affected by the added salt. This may be due to the saturation of 
the ion cloud of the polymer with its counterions and the fact 
that our electric field is imposed externally with no contribution 
from the salt ions. 

To see the effect of the extended electric field and the re- 
lated chain deformation on the capture process, we define the 
monomer approaching probability, MAP(m), of monomer m. 
MAP(m) gives the number of cases (divided by the total num- 




4 5 6 7 

Monomer ID 



Figure 8 Monomer approaching probability (MAP) for a chain length 
of N = 10 for three selected values of the voltage. The inset shows in 
detail the variation of MAP for the end monomers. Note the ascending 
trend toward greater MAP with increasing voltage. Error bars give the 
standard deviation of independent simulations. 



ber of polymer approaches) in which the monomer m is the 
first monomer of the chain which reaches a distance of the 
pore smaller than 2.5(7, starting from an equilibrated config- 
uration. The obtained MAP is plotted in Fig. 7a versus the 
normalized monomer ID, m/N. The plot clearly shows that the 
end monomers have a higher probability to reach the pore en- 
trance as compared to other monomers of the chain. We remark 
that this effect is significantly more pronounced than the well- 
known wall-induced enrichment of the end monomers in the 
vicinity of a planar wall 50 . We thus attribute the main contri- 
bution to this behavior to the deformation of the chain in the ex- 
tended electric field. One can also see that the relative behavior 
of the MAP is nearly the same for all the chain lengths inves- 
tigated. When we compare this result with the rather smooth 
variation of the distribution of the end monomer with distance 
from the polymer's center of mass in its relaxed configura- 
tion 51 , we find a significant deviation which mainly arises from 
the interaction of the electric field and the charged polyelec- 
trolyte. The general behavior of the MAP seen in Fig. 7a is 
consistent with experimental reports on the enhanced capturing 
probability of end monomers 31 ' 32 . 

Next we investigate the effect of counterion distribution. For 
this purpose, we remove the explicit counterions from the simu- 
lation but attribute a constant effective charge (which is length- 
dependent in the investigated N -range) to all the monomers. 
This effective charge is estimated from equilibration simula- 
tions with explicit ions, using the method introduced in a pre- 
vious study on poly electrolytes 46 . 

In Fig. 7b, the monomer approaching probability is plotted 
versus the normalized monomer ID. One can see that elimi- 
nating the freedom of counterions to move and rearrange in 
the inhomogeneous electric field decreases the MAP of the end 
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Figure 9 Distribution of condensed counterions along the chain's 
backbone for different chain lengths, N, for the case of V* = 10 and 
no added salt. The monomer ID (ra) is normalized by the chain length. 

monomers especially for longer chains. One possible explana- 
tion for this observation is that the monomer-counterion bond- 
ing strength depends on the position of the monomer along 
the chain's backbone. While the counterions close to middle 
monomers are more confined and thus cannot detach so easily 
from the chain, it is relatively easy for the counterions close 
to the end monomers to leave the chain. This increases the 
effective charge of the end monomers, thereby enhancing the 
alignment of the chain along the electric field lines. This inter- 
pretation is in line with a previous report on the distribution of 
the condensed monovalent counterions on the chain 52 . 

As seen in Fig. 8 and the corresponding inset, increasing the 
voltage enhances the MAP for the end monomers, while it de- 
creases the MAP for the inner monomers. Again, this is in line 
with the interpretation that the electric field tends to stretch the 
polymer. The larger the applied voltage, the larger this stretch- 
ing and alignment effect, which increases the MAP of the end 
monomers but decreases that of the monomers further away 
from the both ends. This is also in agreement with those studies 
which report a bias for unfolded translocations that increases 
with the applied voltage 14 ' 32 . 

In order to test this idea, we have counted the number of 
condensed counterions along the chain, when it is aligned with 
the electric field lines. We assign the ID number to the end 
monomer which is closer to the pore mouth. As it is visible 
from Fig. 9, the symmetry in the number of condensed counte- 
rions along the chain is broken and there are more condensed 
counterions around the monomers away from the pore. 

In addition to the above discussed effects, long range hy- 
drodynamic interactions can give rise to collective motion and 
thereby influence the translocation process. Focusing on this 
aspect, we have also performed a series of simulations without 
HI so that a comparison with the results presented so far may 
allow to identify the effects arising from HI. In these new set 



Figure 10 Effect of hydrodynamic interactions on polymer deforma- 
tion. The radial extension of the chain is depicted versus radial dis- 
tance from the pore for two different chain lengths of N — 20 and 40. 
The labels WH (NH) stand for simulation results obtained with (with- 
out) HI. V* = 10 and there is no added salt. The right panels show the 
MAP for one half of the chain. 



of simulations, we substitute the lattice Boltzmann solver by 
a Langevin thermostat and thus eliminate the long-range hy- 
drodynamic interactions. We set the friction parameter of the 
thermostat to Fq = 11.8 in order to match the single-particle 
mobility of the Langevin system to the one with full HI. Fig- 
ure 10 compares the results obtained from these simulations 
to those including hydrodynamic interactions. As seen from 
this figure, the effect of hydrodynamic interactions on chain 
deformation is quite small. However, a survey of the monomer 
approaching probability (right panels of Fig. 10) shows that ne- 
glecting HI leads to a significant increase of the MAP of end 
monomers, reaching values as high as 0.48. The fact that such 
a high MAP for end monomers is not in agreement with exper- 
iments 32 shows that HI are quite important for a reliable study 
of translocation dynamics, even though its effects on chain de- 
formation may be quite hard to resolve. 

So far we have only considered inhomogeneous electric 
fields. It has, however, been reported that a polyelectrolyte 
with counterions in a homogeneous electric field greater than 
a critical value, E c , deforms and elongates in the direction of 
the field 52-54 . Now, sufficiently close to the pore, the inho- 
mogeneous electric field may become larger than E c so that, 
at least in principle, the field strength itself, ignoring its in- 
homogeneous nature, can be the main reason for the observed 
deformation. It is, therefore, interesting to study the effects on 
the translocation process which arise from the inhomogeneous 
character of the electric field. 

Since the previous set-up is not suitable to study the effect 
of all the above factors (electric field and its gradient, coun- 
terion rearrangement and HI), we devise a simulation setup 
which allows us to study the effects of both a uniform and a 
non-uniform field. For this purpose, we build a periodic box 
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of size L x x L y x Lz = 103 x 69 x 69 without a pore and put a 
polymer of length 7Y = 30 in the middle of the box. After 10 6 
equilibrium steps, a homogeneous electric field of E = —Eqx 
is switched on and the polyelectrolyte starts to move in the 
positive x direction. This homogeneous field is replaced by a 
spatially varying field after the chain's center of mass has trav- 
elled a distance equal to d = 270 (note that this includes several 
passages through the periodic cell). The inhomogeneous fields 
then acts until the end of the simulation, during which d reaches 
values as large as 340. We estimate Eq using the previous sim- 
ulations as average of the converging field along the centerline 
of the pore from the capture radius to the entrance. This yields 
typical values around E$=E w 0.5. Other parameters and con- 
ditions of these simulations are the same as the previous ones. 

Simulations are performed for four different situations: 40 
simulations without HI and counterions (NH-NC), 40 simula- 
tions without HI but with counterions (NH-WC), 120 simula- 
tions with HI but without counterions (WH-NC) and finally 100 
simulations both with HI and counterions (WH-WC). 

Results obtained from these simulations are shown in 
Fig. 11, where the extension of a polymer (N — 30) along field 
lines is plotted versus the position of the polymer's center of 
mass. Interestingly, regardless of the presence or absence of 
hydrodynamic interactions, a uniform electric field can deform 
a charged polymer only if counterions are present in the sys- 
tem. This best seen by considering the data for d < 270, where 
the applied field is homogeneous. An inhomogeneous field, 
on the other hand (d > 270), can deform a charged polymer 
both in the presence as well as in the absence of surrounding 
counterions. In this context, hydrodynamic interactions seem 
to play a minor role as they only influence the magnitude of 
the effect slightly. Interestingly, the data with HI exhibit sig- 
nificantly larger fluctuations despite the fact that the number of 
simulations in this case is larger than those without HI. We at- 
tribute this effect to the Hi-induced anomalies described in lit- 
erature 54 ' 55 . These anomalies also show themselves in a small 
elongation perpendicular to the field lines of the polymer in the 
absence of counterions (negative slope for d < 270 in Fig. 1 lb). 

Moreover, as the bare polyelectrolyte (no counterions) 
reaches the converging field, it experiences an elongation which 
is significantly larger in the absence of HI as compared to 
the case with hydrodynamic interactions (compare c and d in 
Fig. 11). This shows that, in the absence of counterions, field 
gradient effects are considerably overemphasized if one does 
not take account of hydrodynamic interactions. On the other 
hand, in the more realistic case, where counterions are present 
in the system, the absence of HI only leads to a slight increase 
of the deformation behavior both in a uniform as well as non- 
uniform electric field. 
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Figure 11 Extension of 
a polymer (N = 30) 
along field lines versus 
the position of the poly- 
mer's center of mass, 
(a) In the absence of 
HI (NH), polymer elon- 
gates in a homogeneous 
field (d < 270) only if 
counterions are present 
(WC). In contrast, an in- 
homogeneous field (d > 
270) is able to deform 
the chain even in the 
absence of counterions 
(NC). (b) Similar data as 
in (a) but with HI (WH). 
The general behavior is 
essentially the same as in 
(a). 



4 Conclusions 

In this paper, we study DNA deformation before reaching the 
pore in translocation phenomenon and investigate its effect on 
the capture process by taking into account the electrostatic and 
hydrodynamic interactions. An analytical form for the electric 
field near the pore is obtained by approximating the pore shape 
with a one-sheeted hyperboloid. An approximate expression 
for the electric potential, as a function of the applied voltage 
and the pore characteristics, is also given for the entire whole 
space. This provides us with sufficient information to find a 
theoretical description for DNA elongation. By applying the 
Long's electrohydrodynamic equivalence and Randall's kine- 
matic model for deformation, the velocity gradient tensor and 
hence the extension/compression axes of DNA during the cap- 
ture process are determined. It is then shown that, before en- 
tering the pore, the molecule slightly stretches along the sym- 
metry axis of the pore and, after passing through the pore, the 
stretched DNA experiences a compression along the same axis. 

Using a coarse-grained bead-spring model and considering 
the hydrodynamic and electrostatic interactions, we then study 
via a hybrid lattice Boltzmann-molecular dynamics (LB-MD) 
approach, the deformation of a DNA molecule in the vicinity of 
the pore and discuss its relation to the dynamics of transloca- 
tion. The extension-strain curves show the effective DNA de- 
formation before reaching the pore. It is argued that this helps 
the end monomers to find the pore entrance more easily, thus fa- 
cilitating the translocation process. Our simulations show that 
nearly always, DNA approaches the pore with one of its ends 
closer to the pore and hence most of DNA approaches result in 
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successful translocations. 

Furthermore, by a comparison of simulations with explicit 
counterions with those where counterions are replaced by an 
effective monomer charge, we provide evidence that such a 
mean-field model for the chain charge significantly underesti- 
mates the approaching probability of the end monomers. These 
findings are shown to be in line with experimental studies of 
the capture process. 

One important issue addressed in this study regards the in- 
terplay between hydrodynamic interactions, presence of coun- 
terions and electric field gradients. A result of these studies is 
that, in the absence of counterions, a strong coupling between 
hydrodynamic interactions and electric field gradient occurs. 
In particular, field gradient effects are overemphasized if both 
counterions and HI are absent. On the other hand, if one prop- 
erly takes account of counterions, hydrodynamic interactions 
lead to a more accurate description of the system behavior but 
do not significantly change the overall picture. 

Finally, it is worth mentioning that our proposed electric field 
makes a connection between two distinct fields of research in 
nanopore translocation: (1) Studies in which the access resis- 
tance of a nanopore is investigated and the electric potential 
is assumed to drop mainly inside the pore 30 ' 34 and (2) Those 
in which the effect of the electric field outside the pore on the 
polyelectrolyte diffusion is studied and nearly always the spa- 
tial continuity of the electric field or potential at the pore mouth 
is neglected 12 ' 56 . 
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